Relativistic many-body calculation of low-energy dielectronic resonances in Be-like 

carbon 
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We apply relativistic configuration-interaction method coupled with many-body perturbation the- 
ory (CI-fMBPT) to describe low-energy dielectronic recombination. We combine the CI-(-MBPT 
approach with the complex rotation method (CRM) and compute the dielectronic recombination 
spectrum for Li-like carbon recombining into Be-like carbon. We demonstrate the utility and eval- 
uate the accuracy of this newly-developed CI+MBPT+CRM approach by comparing our results 
with the results of the previous high-precision study of the CIII system [Mannervik et al., Phys. 
Rev. Lett. 81, 313 (1998)]. 

PACS numbers: 34.80.Lx,31.15.Ar,31.25.-v 



I. INTRODUCTION 

One of the important atomic processes governing ionic 
charge abundances in plasmas is dielectronic recombina- 
tion (DR). The DR process is a two-stage reaction with 
a formation of an intermediate doubly excited ion and a 
subsequent relaxation via photon emission, 



^94 



^(9-1) + 



photon. (1) 



Due to the importance of DR in plasma processes, there 
has been a large body of systematic experimental and 
theoretical work on DR. The present status of the field 
is reviewed in Ref. 

An excellent review of current theoretical methods of 
treating DR may be found in Refs. [H, 0|- DR calcu- 
lations were carried out using configuration-interaction 
(CI), multi-configuration Hartree-Fock (MCHF), and 
techniques of many-body-perturbation theory (MBPT) 
with the help of standard codes, such as GRASP ^Sj, 
CIVS [|, MCHF code g, AUTOSTRUCTURE Q, and 
others. For electron temperatures > 100 eV, there 
is a good agreement between calculated and measured 
DR rates. However, for Tg < 10 eV there are significant 
disagreements between theory and experiment (see, e.g, 
Refs. [ll,0,ll|)- These discrepancies are usually attributed 
to theoretical inaccuracies in the positions of low-energy 
resonances Er < 1 eV. Even small shifts of such reso- 
nances to the lower energies lead to underestimating the 
DR rate [|. 

The DR process is a resonant process: cross-section 
spikes at electron kinetic energies that are resonant with 
internal transitions between bound ionic states. As a 
result, the DR rate coefficients, entering, e.g., plasma 
ionization stage calculations, are exponentially sensitive 
to uncertainties in energies of resonances, Er- Because 
of this exponential sensitivity, there is an outstanding 
and practically-relevant problem: a reliable description 
of the DR at low temperatures. This problem has been 



highlighted, for example, by Savin et al. These au- 
thors write, "the single greatest challenge facing modern 
DR theory is accurately calculating the resonance struc- 
ture for the low collision energies needed to calculate low- 
temperature DR. " prl ] Compared to high energies (where 
a simplified Rydberg-like description suffices) , at low ex- 
citation energies the positions of involved atomic reso- 
nances become sensitive to many-body correlations. Solv- 
ing the correlation problem accurately is a challenging 
task, and the existing approaches have difficulties in re- 
liably describing the low-temperature DR. 

The most accurate method to date in treating the low- 
temperature DR is the relativistic many- body theory by 
the Stockholm group (see, e.g., [lOl - flSj and references 
therein). Our present approach shares essential elements 
with this highly-successful method: although our compu- 
tational toolbox is different, it is also based on the many- 
body theoretical treatment and it is ab initio relativistic. 
There is, however, an important difference: all the enu- 
merated calculations by the Stockholm group have final 
ions, y4.'^'?~^)+ with two electrons outside a closed-shell 
core. Our methodology is more general and allows one 
to investigate systems with a larger number of valence 
electrons. The goal of the present paper is to evaluate 
the utility and the accuracy of our approach by compar- 
ing our results with the benchmark data of Ref. flT'l . 

The paper is organized as follows. In Section |TT] 
we present a discussion of basic formulas of DR (Sec- 
tion [H^, the CI+MBPT approach and the complex ro- 
tation method. Specifics of the calculations and numeri- 
cal results for Be- like carbon are presented in Section Hill 



II. METHOD 

A. Dielectronic recombination 

We start by formalizing the DR reaction, Eq. ([IJ, 
and recapitulating well-known results for the DR cross- 
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section and rate coefficient. In an independent parti- 
cle picture, tlie incident electron of energy e excites one 
of the bound electrons of the ion {uala nth) 
and at the same time the initially free electron is cap- 
tured into one of the excited orbitals n£ of the target 
ion: this forms a doubly-excited ion [A'^^^^-'+j . This 
intermediate step is concluded by a radiative decay to 
a final state below the ionization threshold. The theo- 
retical description requires two distinct ingredients: cap- 
ture (auto-ionizing) amplitudes due to electron-electron 
interactions and transition amplitudes due to the pho- 
ton bath. In the isolated-resonance approximation (valid 
for the commonly encountered situation when there are 
no overlapping resonances of the same symmetry), the 
DR cross-section due to an individual resonance may be 
parameterized by the Lorentzian (see, e.g., [T3| ) 
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(e) S'r X — 



r./2 



(AEr - ey + (r J2)' 



(2) 



Here r labels the specific resonant state, AEr is the en- 
ergy of the resonance with respect to the initial state 
\i) of the ion, and is the width of the resonant 
metastable state. For a given range of incident elec- 
tron energies, the total cross-section is obtained by sum- 
ming over the resonances falling within that range. The 
strength, Sr, is given by 



9r 



A°- 4'' 



2rae AEr g, + Aj!" 



(3) 



Here gr and gi are the multiplicities of the resonant and 
the initial ionic states. The radiative decay rate from the 
resonant doubly-excited state, A™"^, is the conventional 
Einstein A coefficient for spontaneous emission summed 
over all possible final recombined states. Usually, limiting 
the radiative decay channels to the electric-dipole tran- 
sitions is sufficient. Finally, A°;^j. is the capture (inverse 
of the Auger process) rate to the resonant doubly-excited 
state and A'^ is the total autoionization rate. The au- 
toionizing and radiative rates account for the total width 
of the resonance, F^ = A™'' -I- A^. 

The above formulation is an approximation: it omits 
radiative recombination (RR), possible effects of inter- 
ference between the RR and DR amplitudes, interfer- 
ence between nearby resonances, modification of the 
Lorentzian profile in the vicinity of the threshold, etc. 
While being approximate, this treatment, however, is 
known to be of a sufficient quality for practical calcu- 
lations :2]. 



beyond the mean- field, Dirac-Hartree-Fock (DHF), po- 
tential) Coulomb interaction between the electrons is 
treated as a perturbation and one applies machinery sim- 
ilar to the textbook stationary perturbation theory. The 
MBPT, when accompanied by a technique of summing up 
important chains of higher-order diagrams to all orders, 
produces excellent results for atoms and ions which have 
only one electron outside closed shells (see, e.g. p^[l6|V 
Even for an atom as heavy as neutral Cs (55 electrons), 
the modern ab initio many-body relativistic techniques 
demonstrate an accuracy of ~ 0.1% for removal energies, 
hyperfine structure constants, and lifetimes [l3]- 

For atoms with two or more valence electrons the 
MBPT has significant difficulties. The dense spectrum 
of many-valence-electron atoms leads to small energy 
denominators and the perturbative treatment of the 
valence- valence correlations breaks down. An adequate 
technique in this case is the CI method. Combining CI 
and MBPT allows one to treat both valence- valence and 
core-valence correlations to high accuracy. 

The CI-I-MBPT technique has been described in detail 
in our previous works |18l-(22{ . Below we briefiy reiterate 
its main features. For simplicity, let us consider a sys- 
tem with two valence electrons outside a closed shell core 
(e.g., a Be-like C'^"'" ion). The calculations start from the 
V^~'^ approximation [21], which means that the initial 
DHF procedure is carried out for the closed-shell ion, 
with the two valence electrons removed. The effective 
CI Hamiltonian for the divalent ion is the sum of two 
single-electron Hamiltonians plus an operator represent- 
ing interaction between valence electrons: 



i/^" =/ii(l) + /ii(2) -1-/12(1, 2). 



(4) 



The single-electron Hamiltonian for a valence electron 
has the form 

hi^ha + t,i, (5) 
where ho is the relativistic DHF Hamiltonian: 



ho — ca ■ p + — l)mc 



V 



N-2 



(6) 



and Si is the correlation potential operator describing 
correlation interaction of the valence electron with the 
core [18]. 

Interaction between valence electrons is the sum of the 
Coulomb interaction and two-particle correlation correc- 
tion operator E2: 



B. CI+MBPT method 



/l2 = 



|ri -r2| 



^2(1,2), 



(7) 



Our calculations are performed using a method which 
combines the configuration interaction (CI) technique 
with many-body perturbation theory (MBPT). MBPT 
provides a systematic prescription for solving the atomic 
many-body problem [14| . Basically, the residual (i.e.. 



Qualitatively, S2 represents the screening of the Coulomb 
interaction between the valence electrons by the core elec- 
trons [ll|. 

A two-electron wave function VE* for the valence elec- 
trons has a form of expansion over single-determinant 
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wave functions 

vl/ = ^c,$,(l,2), (8) 

i 

where $i are the Slater determinants constructed from 
single-electron valence basis states calculated in the 
V^~~'^ potential. The coefficients as well as two- 
electron energies are found by solving the matrix eigen- 
value problem 

(iJ^* - E)X = 0, (9) 

where H°f = {^,\H''^\^j) and X = {ci, C2, . . . , c„}. 

To calculate the correlation correction operators Si 
and S2 we use the second-order MBPT. 

Technically, one needs a complete set of single-electron 
states to calculate S and to construct two-electron basis 
states for the CI calculations. To this end, we gener- 
ate a finite basis set using B-splines and the dual-kinetic 
balance method ^ for the V^'"^ DHF potential. 

By diagonalizing the effective Hamiltonian we find 
the wave functions; they are further used for comput- 
ing atomic properties such as electric-dipole transition 
amplitudes. To compute matrix elements we apply the 
technique of effective all-order ("dressed") operators. In 
particular, we employ the random-phase approximation 
(RPA) . The RPA sequence of diagrams describes a shield- 
ing of the externally applied field by the core electrons. 

The CI-I-MBPT calculations for a system with n va- 
lence electrons follow the very same scheme, with an ex- 
pansion over Slater determinants for n electrons. Again 
the strong interaction between the valence electrons is 
treated within CI, while the core- valence interactions are 
taken into account in the MBPT framework (intermedi- 
ate states in the operator E include core excitations). 

In principle, the effective Hamiltonian H'^^ for sys- 
tems with n > 2 valence electrons includes three-particle 
operator S3, whose computation is very costly. In 
Refs. [1^, [20] this operator was calculated for neutral 
Thallium and respective contribution was found to be 
negligible. One can expect that this conclusion should 
hold for all systems with three-four valence electrons. 
For combinatorial reasons the relative role of the three- 
particle interaction E3 rapidly grows with n and one may 
need to include it for systems with n > 4. 

C. Complex rotation method 

Using the described CI-I-MBPT method we can find 
spectra of multi-valent ions. This section connects the 
CI-I-MBPT method to the DR problem. A straight- 
forward computation within the CI-I-MBPT has certain 
problems, discussed below. These problems, however, 
are elegantly solved using the complex rotation method 
(CRM): a relatively minor modification of the CI-I-MBPT 
method allows us to compute positions and widths of the 
dielectronic resonances. 



What are the difficulties? 
(i) The CI-hMBPT method starts from a finite set of 
single-particle states (orbitals) computed in a spherical 
cavity of radius R. The entire continuum in practice is 
approximated by 20-30 orbitals per partial wave; their in- 
dividual energies depend on R. The DR doubly-excited 
resonance states are embedded into the continuum. In 
some cases, depending on i?, the resonance state may be- 
come degenerate with the quasi-continuum states. This 
leads to a requirement that the model CI space includes 
both the bound and the quasi-continuum many-particle 
states (otherwise the perturbative treatment may break 
down due to small energy denominators) . How would we 
separate the doubly- excited resonance states of interest 
from the background quasi-continuum? 

(ii) Straightforward computation of the capture 
(Auger) rates starting from the Fermi golden rule re- 
quires continuum wave-functions. Because we start with 
the box quantization, we cannot easily generate the con- 
tinuum orbitals of a prescribed energy (in principle, this 
is possible by tuning the radius of the cavity, but this is 
not a very practical solution) . How do we determine the 
autoionizing rates without knowing the scattering states? 

Both difficulties are elegantly solved within the CRM 
framework. 

The complex-rotation method is well established and 
has been employed in atomic physics and quantum chem- 
istry for several decades (see, e.g., a review [23| and ref- 
erences therein). Previously, the CRM was successfully 
applied to the DR problem by E. Lindroth and collabora- 
tors In the CRM, the radial coordinate is scaled 
by a complex factor e*^ 

r re'^ , (10) 

9 being an adjustable parameter. For example, the radial 
Dirac equation becomes 

{V (e'^r) + c2) P„„ {r) + 

Ce"'" ( T" ~ - ) Qnn (r) = SnKP-an (?") , 

\ar r J 

-ce-^' {^ + -]Pn. {r) + 
\dr r J 

{V {e''r) - C^) Q„, (r) - SnnQnn (r) . 

Here P and Q are the conventional large and small ra- 
dial components of the Dirac bi-spinors. For a point- 
like nucleus V (e'^r) — —e^^^Z/r. For the DHF poten- 
tial, the dependence is more complicated as it involves 
^-dependent core orbitals and requires a self-consistent 
solution. 

We have implemented the CRM method for the DHF 
equation using the finite-basis set technique. Technically, 
we employed an expansion over B-splines and the dual- 
kinetic-balance method [2^ to avoid the so-called "spu- 
rious" states. Representative numerical results for the 
Si/2 symmetry are shown in Table HI The generated fi- 
nite basis is suitable for feeding into the CI-hMBPT code. 
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Analytically, the scaling ^TU\\ does not affect energies of 
the bound states, but the eigenvalues of the continuum 
are rotated in the complex plane by —29. Our numerical 
data somewhat deviate from this trend; this is related 
to incompleteness of the spectrum for finite basis sets, a 
result known in the literature. 



TABLE I: Dirac-Hartree-Fock energies for C'^"'' ion in the 
complex-rotation method. Rotation angle 9 = 10°. The finite 
basis set (dual-kinetic balance B-spline basis set) for the 
symmetry consists of 40 orbitals. Cavity radius is i? = 45 
bohr. The £„« < —me(? part of the spectrum is not shown. 



n 






1 


-14.42319872 


-7.09 X 10^" 


2 


-2.365899978 


6.48 X 10"* 


3 


-0.9891728632 


2.61 X 10~^ 


4 


-0.5410495311 


8.00 X 10-^ 


5 


-0.3407231660 


-2.21 X lO-"* 


6 


-0.2339320652 


-9.78 X 10"'' 


21 


519.7385848 


-197.5 


22 


1033.037284 


-380.7 


23 


2013.838233 


-715.3 


24 


3833.246780 


-1296.1 


40 


13970526.56 


-2466750.25 



So far we discussed the one-body problem. It is the 
many-body part of the problem, where the CRM method 
becomes invaluable. When the scaling of the many-body 
Hamiltonian H'^^ is carried out, new, complex, discrete 
eigenvalues of H'^^{9) appear in the lower half of the com- 
plex energy plane. These are the complex values that one 
associates with the resonances: 



E,. — Er, 



(11) 



Here -Eros and r° are the position and the autoionizing 
width of the resonance we are after. For a complete 
set, these complex eigenvalues remain unaffected as is 
varied, while the continuum moves. For a finite basis, 
the Er{9) trajectory in the complex plane "pauses" or 
"kinks" at the physical position of resonances. 

To reiterate, by diagonalizing the complex-symmetric 
CRM-scaled Hamiltonian, we find the 9 trajectories of 
the eigenvalues in the complex energy plane and de- 
duce the positions and widths of the DR resonances. 
Notice that no efforts are needed for computing the 
true scattering states. To efficiently diagonalize the 
complex symmetric matrices characteristic of the CRM 
method, we adopted a specialized Davidson-type eigen- 
solver "JDQZ" 



III. NUMERICAL EXAMPLE: 
BERYLLIUM-LIKE CARBON 

As an ihustration of our CI-fMBPT-l-CRM toolbox, 
we consider DR of the Li-like carbon. The DR pathway 



e" -I- {ls'^2s) 



[C'+ {ls'2p,,nlj)]** 

(...)+ photon. 



Low-energy DR for C^+ was the focus of a combined 
theory-experiment paper [llj ]: we compare our results 
with the results of that work below. 

The calculations were carried out using the relativistic 
basis set with 40 orbitals per partial wave. Representa- 
tive numerical results for the Si/2 symmetry are shown 
in Table HI The correlation operator E was computed 
with partial waves up to Zmax = 5. The CI model space 
was spanned by the 25 lowest-energy virtual orbitals for 
each partial wave up to Zmax = 4. For example, for the 
J = 1, odd symmetry we include all possible 11 angular 
channels. The most important channel is represented by 
the npjn'si/2 and ndju'pj' configurations. The ndjn'pj' 
configurations give rise to leading contributions to the 
low-energy 2p4:d DR resonances. These resonances are 
embedded into the 2sif2£Pj continuum. Including both 
types of configurations allows us to incorporate the in- 
teraction of the DR resonances with the continuum to 
all orders of MBPT and avoid accidental degeneracies. 
While the theoretical formulation of Ref. |ll| is similar 
to ours, their model space excludes the 2snp continuum 
and its effect is taken into account in all-order MBPT. 

Calculations of the CRM trajectories of energy levels 
were carried out for angles 0° — 30° with a step of 1°. 
An example of a trajectory for the 2p4d ^Sq odd-parity 
resonance is shown in Fig. [TJ 

Our numerical results are compiled in Table HH We 
tabulate energies, widths, and strengths of DR reso- 
nances falling within a 0.6 eV range above the ls^2s 
threshold. In this table, we also compare our values with 
the previous theoretical results of Ref. ll|. The ener- 



gies are also compared with the NIST recommended val- 
ues '26*1 . We find that the overall agreement for energy 
positions is excellent and does not deviate by more than 
2-3 nieV. Detailed consideration reveals that the NIST 
recommended values for the position of the 2p4d ^Pj res- 
onance differ both from our and Ref. [ll| predictions 
by as much as 13 meV. Even more strikingly, both our 
and Ref. [Ill predictions disagree with the NIST recom- 
mended value by a very large value of 147 meV for the 
position of the 2pfi4 resonance. 

Among 21 tabulated resonance positions, there is only 
one large disagreement between our work and theory of 
Ref. [ll| : this happens for the 2p4d ^5*0 resonance where 
the two calculations differ by 34 meV. Considering an 
excellent agreement for other 20 resonances between the 
two calculations, such a disagreement may indicate a ty- 
pographical mistake in Ref. jUj. Ref. [11| reports the 
following experimental positions of the resonances: 0.182 
eV for 2pAd'^D, 0.244 eV for the unresolved 2pif^-^F, 
0.438 eV for 2pif^D, and 0.578 eV for 2p4:d^P with ex- 
perimental uncertainty of 5 meV. All of these values are 
in agreement with our theoretical predictions. 
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FIG. 1: (Color online) Illustration of the CI-^MBPT+CRM 
method for locating the 2pAd ^So resonance. The upper (cen- 
tral) panel shows dependence of the real (imaginary) part of 
the energy on the CRM rotation angle 6. The bottom panel 
shows the corresponding trajectory in the complex plane. The 
kink in the trajectory marks the position of the resonance. 



A comparison of resulting autoionizing widths between 
our and Ref. IT] calculations indicates a good agreement 
for broad resonances. The agreement for narrow ( width 
< 0.5 meV) is less satisfying. Experimentally, the width 



of such resonances is determined by the experimental 
convolution function and no definitive conclusions can 
be drawn on the basis of theory-experiment comparison. 
Notice, that the 2pA(PD and 2pAf^^'^F resonances were 
resolved in the experiment llj. The calculated rate in 
Ref. [m was larger than the experimental one by 50%. 
It was not clear whether the source of the disagreement 
was theory or experiment. Perhaps, our disa gree ment for 
the width of narrow resonances with theory may in- 
dicate enhanced sensitivity to details of theoretical treat- 
ment. While there are discrepancies for the autoionizing 
widths of the narrow resonances, our CI-I-MBPT-I-CRM 
strengths of the resonances compare well with calcula- 
tions of Ref. [m (see Table El . 

IV. CONCLUSION 



To summarize, we report developing a new method 
for computing properties of low-energy resonances in 
dielectronic recombination. A high-precision descrip- 
tion of low-energy resonances is particularly challeng- 
ing as it is sensitive to correlations. At the same time, 
uncertainties in the positions of the resonances dras- 
tically affect practically important recombination rates 
in low-temperature plasmas. Our theoretical approach 
is based on combining configuration-interaction method 
with many-body perturbation theory and complex rota- 
tion method (CI-hMBPT-hCRM). The method is ah ini- 
tio relativistic. To gauge the accuracy of the developed 
CI-I-MBPT-I-CRM approach, we computed low-energy 
resonances in Be-like carbon. We find a good agreement 
with the earlier high-precision study by Mannervik et al. 
^11]. While here we studied a divalent ion, our developed 
methodology and computational toolbox is well suited 
for exploring resonances in systems with several valence 
electrons outside a closed-shell core. 
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